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I. INTRODUCTION 



Phase transitions associated with the spontaneous breaking of symmetries in gauge the- 
ories are expected to have played a crucial role in the history of the early universe. In 
particular they imply the production of topologically stable defects Such defects, pro- 
duced at the energy scale of grand unification, lead to a panoply of cosmological effects. 
The most important consequence of their existence, however, is their potential to create en- 
ergy distribution anisotropics, necessary for seeding the large scale structure of the present 
universe, and compatible with those observed in the cosmic microwave background 0. 

Different classes of defects lead to distinct cosmological consequences. The most dra- 
matic is that walls and monopoles are known to come to dominate the energy density of the 
universe, if inflation does not occur subsequently to their formation or some other mech- 
anism intervenes to enhance their annihilation. Strings by contrast, under certain general 
conditions, seem to be entirely viable on their own. The fundamental difference arises from 
the fact that a network of cosmic strings possesses natural mechanisms to reduce its own 
contribution to the total energy density. Moreover, this is done in a way such that the 
effective evolution of the energy fraction in strings decreases in time faster than would be 
expected from the expansion alone and scenarios of string domination are thereby naturally 
precluded. 

The mechanisms that allow for the viable evolution of a string network are essentially 
motivated on topological grounds and are consequently expected to be model independent 
lUH]. They consist of a two step process. Firstly, when two segments of a long string intersect 
a closed loop of string will be formed. This happens when the two segments exchange ends 
after collision, which in turn is a consequence of winding number (i.e., vortex topological 
charge) conservation on the plane. This is designated intercommuting. Finally, because 
closed loops of string are not globally topologically stable they can radiate away the energy 
trapped in their field configuration as gravitational waves and shrink until they disappear. 
The final stage of collapse is probably characterized by vortex annihilation and is expected 
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to produce extremely energetic cosmic rays. 

The realization of this scenario rests crucially upon the efficiency of the intercommuting 
and subsequent string separation. Strings, however are know to experience interactions 
which depend, in general both on the underlying field theory and on the specific region of 
its parameter space. In its simplest and most widely used model strings are the classical 
non-trivial solutions of the Abelian Higgs model. In particular, for certain ratios of coupling 
constants, they are analogous to the vortex solution in type II superconductors and to 
vortices in superfluid 4 He. This latter case corresponds to the vanishing of the gauge sector 
in the model and the corresponding string solutions are naturally known as global strings. 

Typically, when estimating how string networks evolve, numerical simulations only invoke 
the underlying field theory in the initial string formation. Thereafter, they are taken to obey 
the equations of classical Nambu-Goto strings, interaction-free but for empirical assumptions 
about intercommuting. For type II and global strings these assumptions are buttressed by a 
number of numerical studies. They reveal, in particular, that strings indeed exchange ends 
and separate in these cases. Studies regarding the outcome of collisions of type I strings, 
which feel an attractive interaction regardless of orientation, are much scarcer ||. On the 
other hand, analytical studies concerning the derivation and generalization of such range 
of behavior having the field theory as starting point are scarce and tend to concentrate on 
specific regions of parameter space f||7j]. In this paper we attempt a more general analytic 
analysis. 

The article is organized in two parts. The longer is concerned with classical vortex 
interactions. The first sections essentially review the approximate field solutions for vortices 
and derive the interaction potential between two vortices or strings, for all values of coupling 
constant ratios and arbitrary orientations. Models for string production typically lead to 
initial high string densities. The latter part of the paper examines the initial stability of 
string networks, both classically and quantum mechanically. This part of the paper is rather 
more speculative. Type I and type II strings are a consequence of theories that, essentially, 
undergo either first-order or second-order transitions at the time of string formation. We 
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make a preliminary attempt to relate this to the nature of string forces. 



II. THE INTERACTION ENERGY BETWEEN TWO STRINGS: 



GENERALITIES 



Prior to considering multiple string solutions we briefly review the isolated string field 
solutions, with particular attention to the field behavior at large distances from their axes 
of symmetry. Although well-known we need the details for later sections. 

£7(1) strings are the extrapolation to three dimensions of the well-known vortex solution 
of the Abelian Higgs model ||. These are finite energy non-trivial static solutions on the 
plane. The Lagrangian density takes the form 



where e is the gauge coupling constant to the scalar field, A the scalar field self-coupling 
constant and r) the vacuum expectation value of the modulus of the scalar field in the broken 
phase, as usual. 

The finiteness of the energy alone determines the asymptotic form of the fields at infinity. 
In order for the scalar field potential to vanish at infinity, 4> is constrained to lie on a circle of 
radius 77. Its magnitude becomes then well defined but its phase remains arbitrary. On the 
other hand, the modulus of the covariant derivative must also vanish asymptotically which 
forces the gauge field to assume the form §| 



This also ensures that the energy contribution arising from the field strength tensor will be 
finite. Moreover, if we work in polar coordinates, we see that the problem of parameterizing 
the behavior of the scalar field at infinity reduces to the mapping of a circle, in coordinate 
space, onto another, in field space. These mappings, as is well known, fall into an infinite 
number of homotopy classes, each labeled by an integer n G Z. Each solution, characterized 



£ = + \\{d, + ieA,)ct>\ 2 - ^(|0| 2 - V 2 ) 2 



(1) 



lim\g\^ +00 A e (x) = — di(f>(x). 



(2) 
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by a given n, cannot in turn, be continuously deformed into another with m 7^ n, without 
changing the boundary conditions for the fields at infinity. Since this would imply expending 
an infinite amount of energy strings are topologically stable. 

Even though the above arguments supply us with with some information about vortex 
solutions the way in which the fields tend to their asymptotic forms is only revealed by 
studying the Euler-Lagrange equations. These can be solved numerically in all detail but, for 
certain purposes, this does not suffice. Such is the case when we want to discuss analytically 
the detailed energetics of a multiple string field configuration. 

Exact single solutions exist only for the case of critical coupling (b — j- — 1) and positive 



winding number ||10|| . This is a rather special case since the interactions between vortices 
disappear and the equations of motion for the fields become first order; the well-known 
Bogomol'nyi equations. In general, however, approximate solutions in certain regimes can 
be derived. In order to show this we take the Euler-Lagrange equations, assuming cylindrical 
symmetry, in the radial temporal gauge A r = A Q = 



v r V- 



rj '2- ■> ^ = ° ( 3 ) 
VxVx4 s + e^ 2 (eA B =0, (4) 



where ip is the modulus of the scalar field and n its winding number. The field <fi is therefore 
assumed to be of the form </>(r, 9) = Lp(r)e md , where 9 is the angular co-ordinate on the 
plane. To probe how the fields approach their asymptotic values at infinity, one can expand 
ip around 77 as <p(r) = rj — /(r), where / is an auxiliary field. Then the Euler-Lagrange 
equations become 



V 2 / = m 2 s f + f 



e 2 Q 2 ~ \f (37? " /) 



e 2 vQ 2 (5) 



V x (V x Q) = m 2 A Q - e 2 (2 V f - f 2 )Q } (6) 

where m 2 s = rj 2 X, = r/e and Q = Ag — — . 

Equation (6) has an obvious solution for small /, when the second term, on the right 
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hand side, can be neglected in face of the firstQ 



M r ) = k A Kx{m A r), (7) 

er 

where k A is a constant and Kx{m A r) the modified Bessel function of order 1. Equation (7) 
shows us how Ag approaches its asymptotic value. This result allows us to estimate the 
corresponding behavior for the scalar field. Neglecting the quadratic and cubic terms on the 
fields in (5) we obtain 

f = k s K (m s r), (8) 

where k$ is another constant and Ko(msr) the modified Bessel function of order 0. For large 
arguments the modified Bessel functions of order and 1 have the same leading behavior, 
namely, 

«~&"( 1+0 <^)- < 9) 

The approximate solutions (7) and (8) when taken in the limit (9) indeed reproduce the 



exact asymptotic form of the fields at critical coupling ]T0[ . As we have seen above, however, 
they are valid for all values of the coupling constants in the model. 

For short distances, when (rj — <p) ~ r/, equation (6) becomes essentially that for a free 
gauge field. The corresponding form for Ag is then 

Mr) oc — + 0(r). (10) 

er 

This form is actually compatible with (7), when the small argument for the Bessel func- 
tion is assumed. Moreover, if we impose that the magnetic flux should vanish at the origin 
we will determine k A = nrj ||. This is certainly an acceptable procedure for strongly type 
II vortices. An approximate solution for the scalar field close to the origin can equally be 
found by using (10) in (5). We obtain 



1 This is a special case of the most general solution Ag(r) = — — kAKi(e<pr), when 99 ~ constant. 
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<p(r) oc Jn{^= r ), (11) 

where J n is the Bessel function of order n, which is the usual winding number. This implies 
the well-known small distance behavior for 

^ r)(X l{^ r ) ' (12) 

In particular for n — 1 the field is linear close to the origin. Because the behavior of the scalar 
field in these two quite different limits is not directly relatable, k$ cannot be computed with 
generality. In the rather special case of critical coupling (y = 1), however, the dynamical 
equations for the fields reduce to 

Ti 

d r ip ± (eA e - -) = (13) 
r 

d r ^±^V-r? 2 ) = 0, (14) 

where ± = sign(n). These are the Bogomol'nyi equations for critically coupled vortices 
in polar coordinates. It is easy to see that they allow us to relate / and A e directly and 
consequently also ks and k,A- Using the identity of the two coupling constants we obtain 
ks = \n\r). In principle there is no reason for this result to hold in other regions of parameter 
space other than the good agreement of its consequences with the numerical results, as we 
shall see later. Nevertheless unless otherwise stated we will assume it henceforth. 

We are now in a position to derive the interaction potential between two string seg- 
ments, under the assumption that the solution for the two string field configuration can be 
successfully approximated by the superposition ansatz |TT|, 



(h( \ <P(\r - r x \)(j){\r - r 2 \) 

$(r,ri,r 2 ) = (15) 

Ae(r, ri ,r 2 ) = A e (\r - n|) + A e (\r - r 2 \). (16) 

Here <p an d Ag are the isolated string field configurations of Section 1. In order to 
estimate the interaction energy between the two vortices one then simply substitutes (15) 
and (16) in the energy functional 



E[<j>,A li ] = jdV 



\F, V F^ + l -\d r 4\ 2 + l -\eA e 4> - d e <t>\ 2 + ± 



2 - v 2 ) 2 



(17) 



We then obtain 



E[$, A] = Efa, A,] + E[4> 2 , A 2 ] + J dV 1 - (2V ri x A e { ri ).V r , x .4 e (r 2 )+ 



+ 



if 



(y^) 2 - T?)(V n <p{r 2 )) 2 + (<^(r 2 ) 2 - rj 2 )(V r Mn)) 2 + 



+ 2v9(ri)v?(r 2 )V ri v?(r 1 ).V r2 v?(r 2 ) + 

+ (^(r 2 ) 2 - ^MnfieAein) - ^) 2 + (y^n) 2 - V 2 Mr 2 )\eA e (r 2 ) - ^) 2 + 

ri r 2 



Tli. 



n 2 , 



+ 2y?(r 1 )V(r 2 ) 2 (eA,(r 1 ) - -).{eA e {r 2 ) - -± 



r 2 



+ 



+ 



_(v?(ri) 2 - r/ 2 )(^(r 2 ) 2 - r? 2 ) [^(n) 2 + ^ 2 )(^(r 2 ) 2 + r? 2 ) - 2r/ 4 
Finally, we subtract the contributions due to the isolated vortices. We obtain 
E int [ip,A e ] = J dV 1 - (2V ri x A e (r!).V r2 x A e (r 2 ) + 



(18) 



+ 



^{nf - v 2 )(V r Mr2)) 2 + (^(r 2 ) 2 - v 2 )^rMn)) 2 + 



+ 2ip(r 1 )ip(r 2 )V ri (p(r 1 ).V r2 (f(r 2 ) + 

+ (y?(r 2 ) 2 - 7 7 2 )^(r 1 ) 2 (eA,(r 1 ) - ^) 2 + (y^n) 2 - V 2 Mr 2 ) 2 (eA e (r 2 ) - ^) 2 + 



ri 



r 2 



n l . 



n 2 . 



+ 2v(n) 2 tp(r 2 ) 2 (eA e (r l ) - -).{eA e {r 2 ) - ^) 



ri 



7*2 



+ 



+ 



v 



(V^(ri) 2 - ^ 2 )((^(r 2 ) 2 - r? 2 ) (y^) 2 + r/ 2 )(^(r 2 ) 2 + r? 2 ) - 2rf 



(19) 



This last expression can also be written in terms of the auxiliary field /. Assuming 
cylindrical symmetry the energy per unit length of the strings becomes 



£int = J dS {E x H 2 el x .e^ 2 + m 2 A QiQ 2 ee 1 .ee 2 
- (V r /iV r / 2 epi.e72 + m|/i/ 2 ) 



/ 2 Q 2 + 2/^x^2 + / 2 Q 2 + 2/ 2 Q 1 .Q 2 + 2/ 1 j 2 (Q 2 + Q 2 ) + Af x f 2 Q x .Q 2 



+ A(3/ 2 / 2 + 3/2/0 



»7-g(/i + /2) 



+ gA(/ 1 3 /2 + / 2 3 /i)- 



(20) 



To obtain the detailed behavior for the interaction energy one then simply needs to 



perform the integrations in (19) or (20) assuming given functional forms for the fields. This 
procedure is absolutely straightforward but a few cautionary remarks are in order. 

Firstly, unlike the case of critical coupling, the proof for the existence and uniqueness of 
multi- vortex solutions for general couplings does not exist. In particular, any such solution 
will not be static. However it is known, from experimental evidence in type II superconduc- 
tors that multi-vortex configurations indeed do occur. 

Secondly the above procedure is clearly not exact. The superposition ansatz effectively is 
not an exact solution of the Euler-Lagrange equations (3) and (4). The resulting discrepancy 
on the value of the interaction energy can nevertheless be computed. This can be achieved by 
assuming residual additive terms to (15) and (16) such that the resulting fields would solve 
(3) and (4) exactly. Their behavior and contribution to the energy can then be estimated 



whenever the functional form for the individual vortex solutions is known |12] . In particular 



for large distances away from the strings' axes, when the modified Bessel function behavior 



for / and Q applies, such contribution is of the order ^4- [13]. This is a second order 



subleading effect as we shall verify in the next sections. This whole procedure allows us 
consequently not only to state the validity of the superposition ansatz for a given functional 
form of the fields but also to compute the magnitude of the approximation involved. 

Finally, the integrations in (19) or (20) necessary to obtain the total energy can, in 
general, also only be performed approximately. This will be the subject of the two ensuing 
sections where expressions for the interaction energy are indeed obtained, for all values of 

2 

b = y and arbitrary orientation of the two string elements. 

III. THE INTERACTION ENERGY BETWEEN PARALLEL STRINGS AS A 

FUNCTION OF B 

In this section we derive the integrated expressions for the interaction energy, for all 

2 

values of b = For the sake of simplicity we will restrict ourselves to the case of two 
parallel strings and leave the study of arbitrary geometries to the next section. The problem 
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then reduces to the study of two coplanar vortices. Because different values of b imply quite 
distinct field configurations we will analyze the cases of critical coupling (6 = 1), type II 
(b < 1) and type I vortices (b > 1) separately. 

A. The case b = 1 

In view of the approximate field solutions obtained in Section 1 and the expressions for the 
interaction energy (19) and (20) derived in Section 2, this is the simplest case to approach. In 
general there are two length scales associated with each vortex. They are simply the inverse 
classical masses for the gauge and scalar fields i.e., ta = ra^ and rs = rrig 1 , respectively. 
When the scalar field acquires a non-zero expectation value (< ip >^ r -+oo v) both fields 
become massive and the modified Bessel function solutions for the fields constitute a good 
approximation. This should happen for distances R, measured from the zero of the scalar 
field, such that R> r s . 

These approximate solutions tell us how the fields approach their asymptotic values, 
which in turn correspond to the vacuum of the theory. At distances R larger than r s a 
vortex only perturbs the vacuum by effectively acting as a source for the fields Q and /. 
The interaction between two vortices, separated by a distance d > 2r s then reduces, in the 
first approximation, to the interaction between these two fields. This picture can indeed be 
obtained from (20) if we keep only the terms linear on the fields of each vortex. Then we 
have 



Higher order terms will be much smaller due to the fact that the fields themselves are 
exponentially decreasing. Moreover, their contribution is of the order of the approximation 
involved in adopting the superposition ansatz for the two-string field configuration. 

To integrate (21) we first note that e^.e^, is independent of the point on the plane. 
Then, using the fact that the modified Bessel functions of order zero satisfy 




(21) 
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(V 2 - m 2 ) K (mr) = -2ir6(r), (22) 
in two dimensions, we can write 

E int = lixtf [n^KoirriAd) - \ni\\ri2\K (msd)] + 



+ 



dSnin 2 ^K (m A r 2 )V 2 ri Ko(m A r 1 ) + m A K 1 (m A r 1 )K 1 (m A r 2 )ee 1 .ee 2 ) - 



^K (msr 2 )V 2 ri K (m s r 1 ) + V ri K (m s ri)V r2 Ko(m s r 2 )er 1 .e^ 2 ) . (23) 
Now by noting that 

T\ — dcos(9i) 

v r Ko[mr) = —mKi[mr); \7 ri Ko(mr 2 ) = — m Ki(mr 2 ), 



r 2 

and that 



(24) 



-, -, r 1 -dcos9 1 

e ei .e d2 = e ri .e r2 = , (25) 

we see that the integral in (23) vanishes and the interaction energy becomes exactly 

E int (d) = 2-nrf (mn2K Q (mAd) - \nin 2 \K (m s d)) . (26) 

This shows in particular that, in the limit of our approximations, two vortices at critical 
coupling are free. This is equally true for a pair of parallel cosmic strings, in agreement with 
well-known numerical results 0. 

B. The case b < 1 

Away from critical coupling the two length scales r$ and r A become different and the 
fields are expected to change accordingly. In the specific case when e 2 < A, r A will be larger 
than r$- These two lengths now define a set of three coaxial cylinders, around the axis of 
the string, where the fields may behave in quite a different fashion. Essentially, for distances 
R larger than r A or smaller than r$ the picture of the last subsection remains qualitatively 
unchanged. Indeed, for R < rg the fields must vary substantially from their values at the 
origin to a value not too far from their asymptotic forms. 
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For R > ta, on the other hand, the scalar field should be very close to its asymptotic 
value and consequently will be a slowly varying function of r. Then both fields will have 
acquired their classical masses. For this case the computations of subsection 4.1 are certainly 
valid and it follows that two vortices will repel and a pair vortex-antivortex will attract, with 
magnitude given by (26). 

The behavior in the intermediate region, however, turns out to be the most interesting. 
There R > rs and the fields are expected to acquire mass. This implies, namely, that the 
modified Bessel function form for Q holds. However because itiaR < 1, the Bessel function 
can be approximated by its small argument behavior 

Tl 

Q = —nr)K 1 (m A R) ^ -• (27) 

etc 

When this form is used in the Euler-Lagrange equation for the scalar field, assuming ip = 
rj — f, we obtain 

n 2 1 

^ (m^r) 2 ^ ^(m^r) 4 ^' ^ ^ 

which is the solution for the global string, in the same regime (R > rs) @. 

It should now be clear that, as we take e 2 — > 0, this intermediate region grows to fill all 
space (from R = rs to infinity), and the global string behavior is fully recovered. This is not 
very surprising, of course, since in this limit the gauge field decouples from the scalar field 
and the Euler-Lagrange equation for the latter is that for the theory exhibiting the global 
symmetry alone. What is nevertheless quite interesting is that for an ensemble of type II 
local strings to interact as in the global case they simply need to lie at relative distances 
smaller than r^.Q 

Let us now see in detail how this picture arises from the previous expressions for the 
interaction energy. The part of the energy dependent on / can be integrated to give a 



2 We should keep in mind that this result is a direct consequence of the value taken for !ca, in 
Section 1. 
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contribution which falls with at least m l d3 . Because msd > 1 this contribution will be 
small. On the other hand since the modified Bessel function of order 1 is still a solution for 
the gauge field the computations of subsection 4.1 referring to it will still hold. We then 
obtain, keeping the terms on the gauge field only 

E(d) ~ 2-nrfn 1 n 2 KQ{m A d) ~ — 2n7] 2 n 1 n 2 ln(m A d), (29) 

since ra^i? is small for R < r A - However, our assumptions about the nature of the fields 
leading to (29) were only valid for R > r$, whereas in (26) we implicitly took it to be 
applicable in all space between the origins of each vortex. We must, therefore subtract the 
energy corresponding to R < r s . We then get 

E(d) ~ — 27r7] 2 nin 2 (\n(m A d) — ln(m^rs)) = — 27n] 2 n 1 n2 ln(msd). (30) 

This expression is just what we would expect for the global string, c.f. ||, and ensures the 
consistency of the arguments above. As is well known this results in a strong, infinite rangeQ, 
repulsive (attractive) force for a pair vortex- vortex (vortex-antivortex). 

C. The case b > 1 

In the preceding analysis we were able to compute estimates for the interaction energy 
between vortices under a few quite general assumptions. In particular we learned that the 
proximity of the scalar field to its asymptotic value (traduced in the criterion R > rs) 
introduced a fundamental qualitative change in the behavior of the fields, allowing us to use 
the approximate solutions (7) and (8) with some confidence. The case of type I vortices is 
naturally characterized by the fact that r A < rs- The results for b = 1 for the interaction 
between two vortices should in particular still be valid as long as their scalar cores do 
not overlap, i.e. for (d > 2r$). Then, unlike what happened when b < 1, and because 



In the case of nonvanishing e, m A provides the natural cutoff at large r 
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now itla > rns the scalar field contribution will dominate the energy. Moreover, since 
the corresponding term in the interaction energy (26) has no dependence on the sign of the 
winding numbers this will always result on an attractive force regardless of the actual nature 
of the two vortices involved. 

As may easily be anticipated we will run into a problem when we try to estimate what 
happens in the intermediate region {ta < d < rs)- This can be explicitly shown if we search 
for the behavior of the fields of a single vortex in this region of space. 

In fact, the way we derived the Bessel Function behavior of Q for large distances was to 
assume that tp would be close to r], so that the mass term in (6) would clearly be the dominant 
one. This always holds provided that R ~ r s , and is thus certainly a good approximation 
for the cases above. 

If, on the other hand, R is small enough so that (y? 2 — rf ) ~ — rf, and Q exhibits its 
short distance behavior, then the solution of (3) becomes 

t \ T r msr \ 1 f Tn S r \n foi\ 

f(r) oc Jn{~^) - -^r^ • ( 31 ) 

But this is the behavior we expect for the scalar field very close to the origin, indepen- 
dently of any particular value of the coupling constants. 

However, we also note that as we progress towards R = rs, f approaches rj and the 
Bessel function behavior (7), and (8) start being valid for Q and /, respectively. When this 
happens, Q automatically assumes an exponentially decreasing form since m^i? is large but 
/ should still behave approximately logarithmically since rrisR is small. This should happen 
in a region where the distance to the string's axis is sufficiently close to r s . Consequently 
and unlike what happened for type II strings this approximate behavior for the fields can 
only exist for a thin region of space around a string. Its consequences for the interaction 
energy between the two strings will then only constitute a transient regime between its form 
for large distances and the actual superposition of the two vortices. Nevertheless we can 
compute what the corresponding interaction energy should be. Proceeding in the same way 
as in 4.2, but now only keeping the contribution from the scalar field we obtain 
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E(d) ~ —2iir] 2 \nin2\Ko( y msd) ~ 27r?7 2 1 rz,i ri 2 1 lnfm^d), (32) 

where, again, we implicitly assumed that this regime was present in the whole region between 
the two axes. We should then remove the contribution from the inner region where this does 
not apply. In the absence of another length scale we can only subtract the energy arising 
from R < ta- This however will result in an overestimate (in absolute value) for the final 
result. We then obtain 

E(R) ~ 2irT 1 2 \n 1 n 2 \ln{m A d). (33) 

This always gives rise to an attractive force, as expected. 

At distances smaller than rg, the scalar cores of the two strings will superimpose. This 
results in a field configuration of winding number n\ + n 2 as the distance between the zeros 
of the scalar fields vanishes. Such a transformation of the two string field configuration is 
accompanied by the change in behavior between the modified and unmodified Bessel func- 
tions solution for the scalar field of section 1. This change is dictated by the nonlinearities 
in the corresponding Euler-Lagrange equations which implies the breakdown of the super- 
position ansatz. Close to d — 0, however it is known numerically |J that the interaction 
energy should become approximately constant, signalling the fact that a type I string of 
higher winding number is a stable solution of the Euler-Lagrange equations relative to its 
lower winding number isolated constituents. [Alternatively, type II strings of higher winding 
number are unstable with respect to decay into strings of lower winding number] . 

IV. INTERACTION BETWEEN TWO COSMIC STRINGS: GENERAL 

GEOMETRY 

In the previous section we analyzed how the value of the ratio of the two coupling 
constants in the model b changes the field configurations and derived the interaction energy 
between two vortices, i.e. between a pair of parallel strings, per unit length. 
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For a general network of cosmic strings we know [|T^] that strings wiggle and bend, 
possibly on several scales, and, in general, no two strings are parallel^ We shall expand 
upon this in the next section. 

The assumption of parallel strings, however, conveniently allowed us some comfortable 
simplifications. Because each string is locally (i.e. for any element of infinitesimal length 
along its axis) cylindrically symmetric we could associate the natural frame for the second 
string with that for the first simply by translating it by a distance d, within the same plane. 

This distance was then automatically defined for all pairs of string elements and the 
interaction energy for a piece of string due to the presence of the second is simply that 
for one element times the length. For a general configuration of the two strings both the 
distance and orientation have to be specified for each pair of elements. They will therefore 
weigh differently in the evaluation of the interaction energy between two pieces of string of 
finite length. 

Because of the cylindrical symmetry of the fields composing a string is still locally pre- 
served, to find the interaction energy of one of its elements we need only look for an element 
of a second string intersecting the plane in which it lies. The natural frame for the element 
of the second string will appear in the general case to have its origin at a distance d but also 
to have been rotated so that it no longer lies in the plane of the former. However, because of 
its own local cylindrical symmetry, rotations around its axis of symmetry leave it unchanged. 
This rotation relative to the plane of the first vortex is then generated by two angles only. 
This can be seen in figure 1, where we chose the x-axis along the direction connecting the 
origins of the two vortices. The angles a and 7 parameterize rotations around the x and 
y-ax.es, respectively. 

In order to be able to generalize the results of Section 3 we must be able to understand 



4 Unlike cosmic strings, flux lines in superconductors and vortices in superfiuids are indeed parallel. 
For those cases the results of Section 4 apply. 
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how the change in orientation of the two relative elements affects the interaction energy. 
This will result essentially in a change on the inner products between the unit vectors for 
the directions associated with the natural frame of each string. 

This change can be computed by rotating one of the frames around the x and y-axes of 
figure 1. We now see, in particular, that e^.e^ will simply have the form 

e^-e^ = cos(a) cos(7). (34) 

This allow us to integrate the two terms in (23) to obtain 

E int ( d, a, 7) = 2nr] 2 [nin 2 Ko(m A d) cos(a) cos(7) - \ni\\n 2 \K (msd)] - 

- j dSnin 2 {V ri Ko{m A r 2 )V ri K (m A ri) cos(a) cos(7) - Ki(m A ri)Ki(m A r 2 )e^ 1 .e^ 2 ) 

- KIM (V ri K (m s r 2 )V ri K ('m s r 1 ) - V ri Ko(m s r 1 )V r2 K (msr 2 )er 1 -O] . (35) 

The two residual terms, inside the integral, behave somewhat differently. In the case of 
that arising from the scalar field the change in variable of differentiation on the first exactly 
generates the inner product e^.e^. Because both these terms have opposite signs their sum 
vanishes. 

Such is not the case for the corresponding term in the gauge sector. There we have 

-r~ e zi- e z2 = cos(o;) cos(7) Xi(xi — d) cos 2 (7) + y\ cos 2 (a) 

dr x r x r 2 L 

+ X\Z\ cos(a) cos(7) sin(7) — x\y\ sin(ct) sin(7) cos(7) 

+ y\Z\sin{a) cos 2 (a) , (36) 

which must be compared with 

r cos(a) cos(7)[rri(xi — d) + y\] 



r\r 2 

2 



+ X\Z\ cos (a) sin(7) — x±y± sin(a) cos(a) sin(7) + yiZisin(a) cos(7) . (37) 

In both the above expressions we took a rotation by 7 after a rotation by a. 
These two terms do not cancel in general. Rewriting (35) we finally obtain 
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Eink{d, a, 7) = 2-Krj 1 [nin2K (mAd) cos(a) cos(7) — |ni||n2|-K"o( m s ,( ^)] — 

- [ dSn 1 n2Vr 2 K (m A r 2 )Vr 1 K (mAri) (-^— x 1 (x 1 - d) sin 2 (7) + 
J \rir 2 

+ y\ sin 2 (a) + higher order in a, 7 J cos(a) 003(7). (38) 



We see that the interaction energy is essentially given by the two first terms. In the limit 
of vanishing a and 7 the residual contribution in the integral effectively goes to zero with 
the sine squared of the angles, in agreement with the results of Section 3. 

Expression (38) stresses the quite different nature of the two contributions to the inter- 
action energy. The vector-like character of the gauge field introduces a dependence on the 
relative orientation between the two strings. In particular we see that if one of the strings is 
rotated by 7r the corresponding interaction energy term changes sign. This is equivalent to 
changing the nature of a vortex into an anti-vortex or vice-versa and results from the geo- 
metric nature of the winding number. The term arising from the scalar field is, in contrast, 
insensitive to any particular configuration for the two-string system, as could naturally be 
anticipated. 

It also becomes clear that critical coupling (b = 1), for an arbitrary geometric configu- 
ration, ceases to be the special case when the interactions between strings vanish since the 
angular dependence on a and 7 destroys the balance between the gauge and scalar terms in 
the energy. 

All expressions derived above concern elements of string with infinitesimal length, re- 
sulting from the integration of the interaction energy functional on the plane. To obtain 
the interaction energy for two segments of string of finite length one then has to integrate 
over a series of planes perpendicular to one of the strings. In particular, every time the 
string bends so that any of its elements will lie at an angle larger than ~ it is expected to 
experience self-interactions. This is always the case for closed loops of string which contain 
part of the energy arising from self-interactions, in addition to their tensiorfl. 



5 If the string's tension is defined to be the energy/unit length for an infinite straight string 
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The expressions derived above allow us to compute analytically the interactions for a 
variety of configurations. Let us consider two simple examples. 

For a circular loop of large enough diameter D (D > 2rs as always) it is straightforward 
to compute its self-energy. It is simply 

E self (D) = ^-E(D), (39) 

where E(D) is given by expression (26), taken in the appropriate limit, for a parallel pair of 
string-antistring (|ni| = — 1?7-2 1 ) - 

Another interesting case is to compute the interaction energy due to a circular loop 
of string, in the limit when the distance to a test string segment is much larger than the 
loop's diameter. Then, see figure 2, the interactions are essentially due to a dipole of 
string segments. The corresponding interaction energy is proportional to the sum of the 
Bessel functions coming from the gauge and scalar field contributions. In particular for the 
contribution from the former we will have 

E int oc (K (mr) - K (mr 2 )) , (40) 

where 



r 2 = JD 2 + r 2 + 2Dr cos(0) . (41) 



Now, in the limit of small diameter compared to distance — << 1 and for the strong type 
II small loops of string, when Dm^ cos(0) < 1, we can series expand the Bessel function to 
obtain 



K I mr\ l + ^- + 2— cos(0) I ~ K (mr) 1 + m 2 D 2 cos(0) 



— mD 

The energy can then approximately be written as[| 



cos(#) sin 2 (9) Ki(mr) 

r 



(42) 



6 Here we have neglected the contribution form the scalar field and took the winding number of 
the test string segment to be one. It is also assumed to lie in a plane parallel to that of the loop. 
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E int ~ 2nrj n mD 



cos(0) - -sin 2 (9) 
r 



Kximr) - m 2 D 2 cos 2 (0) K (mr) . (43) 



This expression has two different limits depending on the relative value of m~ A and r. If m^r 
is large then we can take the asymptotic forms for the Bessel functions at large arguments 
(9) to obtain 



E- mt ~ 2nr] 2 J^^e' mAr (m A Dcos(6) - sm 2 (9)— - m 2 A D 2 cos 2 (6)) . (44) 
y 2m at \ r J 

If, on the other hand, r < m~ A ( extremely type II case), we can use the form of the Bessel 
functions for small arguments to obtain 



E int ~ 27TT] 2 



uiaD [ cos(#) — — sin 2 (0)^) — h m 2 A D 2 cos 2 (9) logirriAr) 

r J m A r 



(45) 



We see in particular that in this latter case the leading order term in (45) coincides with the 
usual - potential for a point-like field source in 3 dimensions. A small loop of (quasi)global 
string thus behaves effectively as a monopole when seen at sufficiently large distances. 

V. STRING FORMATION AND STRING FORCES 

In simulations of string production (e.g., see [T4]]) </> field phases are laid down at random 
in a space divided into correlation volumes within which the field phase is held approxi- 
mately constant. In these circumstances strings with winding number \n\ — 1 are produced 
at high density. The laboratory production Jl5| of string defects in superfluid 4 He has sug- 



gested that this model is plausible |16j , at least for the non-relativistic vortices of the global 
U(l) Ginzberg-Landau theory. 

In addition to the role string interactions have in determining the outcome of string 
collisions they may affect the early string evolution if, as expected, strings are produced 
at sufficiently high density. In the remainder of this paper we shall examine the effect of 
interstring forces in a newly created, and approximately static, ensemble of strings. In this 
section we shall investigate whether, with only attractive forces, type I strings can exist at 
high density. This is not a question of how two strings of lower winding number combine 

20 



to form a string with higher winding number, since we expect strings, of whatever type, to 
appear as a tangled mess. However, we shall return to this point later. To see how type I 
strings could exist we adopt a classical toy model possessing some of the characteristics of 
real strings, in which the effect of attractive forces is more transparent. 

Cosmic strings are produced at phase transitions in the G.U.T. era but, nominally in- 
dependently of this, it is known that thermodynamic ensembles of classical string naturally 
display transitions. Most naively, if cosmic strings are treated as non-interacting Nambu- 
Goto strings with modes of arbitrarily high frequency they undergo a Hagedorn transition 
[17]]. In this section we take this classical transition seriously, but we shall be more realistic 



by treating classical strings more as polymers. These again show transitions but, unlike the 
fundamental Nambu-Goto strings, permit the inclusion of string forces at low densities. 

The comments that follow are a straightforward extension of earlier work of ours on 
global strings |H| . Since they are somewhat speculative a detailed recreation of these earlier 



results is inappropriate, and the reader is referred to the literature. 

Strings produced in the way indicated above behave like random walks of high density 
|TJ|. It is well-known that an ensemble of non-interacting random walks at temperature 
T = P' 1 , step length a, can be described by a dual field theory of a free complex field 
2"0fl . The "loop field" xG&) has action (Hamiltonian) pOfl 



S = I dx 

where 



|Vx| 2 + M 2 | X | 2 , (46) 



M 2 cx (47) 

in which e is the string energy/unit length and s = Ola" 1 ) is the entropy/unit length. The 
vanishing of M 2 at e = Ts defines a Hagedorn transition at T = T^, at which string is 
produced copiously. For T > Th the theory is unstable^. 



7 And at which all energy goes into the production of a single string. 
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This is already enough to indicate why, even if the string formation mechanism naturally 
generated strings with higher winding number \n\ > 1, the thermal energy would be sufficient 
to unpeel them to \n\ = 1 strings. The problem is akin to that of adhesion in macromolecules. 
As a first step we adopt the simplifying assumption that, once separated from a string of 
higher winding number, \n\ = 1 strings experience no forces. A straightforward extension 



1 21] of the work of Wiegel |22| shows that, provided the energy required to split an n = 2 
string into two n = 1 strings is less than the cost of creating the string, the Hagedorn 
transition survives at the same temperature T#. Moreover, there is a temperature range 
To < T < Th in which the separation into n = 1 strings is total. We shall not consider the 
problem of higher winding number further. 

The forces between the strings modify (46). Orientable forces are characterized by a dual 
vector field A. It is not possible to mimic our cosmic strings exactly. However, for a strong 
type II theory of static strings at low density, a plausible dual field theory is given by an 
action 

'1 .s 1 



Srr = dx 



I ,VAi) + ^m 2 A A 2 + I (V + if A) X \ 2 + M 2 \ X \ 



(48) 



where rriA — 0, and f 2 = O{0 2 e) is the dimensionless coupling strength. M 2 is as before. 
By the same token, the extension of Su to include a scalar field a, as in the action 
1 



S = dx 



- A (V A A) + -m\A 2 + | (V + if A) X\ 2 + 

+ 1 -(Va) 2 + 1 -m 2 a 2 -g 2 a 2 \x\ 2 + M 2 \ X \ 



(49) 



corresponds to the inclusion of non-orientable sources and hence the possibility of type I 
strings. The coupling g 2 > characterizes the attractive nature of these forces. The simple 
diagonal form ( in terms of a and A) cannot fully reproduce the results of (38). However, 
it does indicate how attractive forces of finite range are permissible at low density. 

To see this, we note that, in the dual theory of (46), (48) or (49), the dimensionful string 
density p — y, the length/unit volume, is given by < |x| 2 > |3 • We see from (49) that, 

2 

because of the attractive forces, the classical cr-field sector is unstable unless \x\ 2 < -^f-, i.e., 
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unless the string density is sufficiently low. We note that long-range wholly-attractive forces 
would always lead to (classical) instability in this simple model. However, the result that 
finite range attractive forces do not destabilize the network at low enough density (large 
enough separation) seems very plausible. Further, insofar as g 2 increases as the strings 
become more strongly type I, the maximum permissible density drops, given the constancy 
of m 2 s . 

To do better would require a greater confidence in (47) that we can expect. However, 
a one-loop approximation to the effective density potential in which the A and a fields are 
integrated over in (49) gives 



V(p) of (50) is sufficient to show why, when arbitrary string orientations are taken into 
account, there is nothing special about the value b = 1. For b near unity, when ms — rriA 
and (presumably) f — g, the effect of the attractive non-orientable third term is to partially, 
but not wholly, cancel the effect of the orientable force term at low density. This is as we 
would expect. 

For strong type I strings the situation is totally unclear, and there is no natural way to 
identify type I and type II strings with first-order (melting) and second-order (Hagedorn) 
transitions, directly. It could be argued that neither should there be, since whether a string 
is type-I or type-II has more to say about the nature of the quantum field phase transition 
than classical transitions. It is to the quantum aspects of string production that we now 
turn. 

VI. STRING FORMATION FROM QUANTUM FLUCTUATIONS 

The main reason for halting our analysis of the dual theory is not so much its incomplete- 
ness as the fact that, initially, cosmic string production arises from quantum fluctuations at 




(50) 
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the phase transition possessed by the underlying local U(l) field theory. Without reference 
to string, in approximate thermal equilibrium, the theory displays a second-order transition 
at the temperature T c = 0(rj) if b << 1. On the other hand, for b >> 1 it is reliably 
expected to display a first-order transition. This accords with our terminology of type II 
and type I strings, respectively. However, there is no reason to believe that the transition 
changes from second to first order exactly at b = 1. 

Whatever the case, in quantum field theory there cannot be a Hagedorn transition, with 
its characteristic maximum temperature. Rather, there is a critical temperature T c above 
which the U(l) symmetry is restored and strings cannot occur. In thermal equilibrium 
at temperature T, this is best seen by working in Euclidean time, in which the fields are 
periodic with period (3 = T _1 . Provided T is larger than the mass scales the 'heavy' (n ^ 0) 
modes in the Fourier series can be integrated out, to give an effective three-dimensional 



theory derived from (1), with action p3|j24 
PS3 = J dx 



~ (V A A) + ~ | (V + ieA) 4>\ 2 - \m\TM + ^ 



+ terms containing Ao(x) + counterterms, (51) 

where A$(x) is the temporal component of the "light" (n = 0) mode component of the 
gauge- A field, and 

-™ 

is the effective scalar mass of the theory. 

There are two points to note about the action (51). Firstly, when A = 0, it is extremised 
by the same (static) vortex solutions as was the action based upon (1), but for the fact that 
the Higgs mass is temperature dependent. The second is that the Higgs potential in (51) 
shows a second-order transition, prior to integrating over the gauge field A, whatever the 
value of b. A first order transition can only occur for b >> 1 as a consequence of gauge-field 
radiative corrections. 

Returning to the first point, we could use the 'classical' solutions to (51) as the basis 
of a dual field theory, as in the previous section. This would now be a dual theory in the 
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more usual sense of the word, a rewriting of the original quantum theory in terms of its 
excitations. The instability of this theory is, in the first instance, characterized again by the 
vanishing of 

U = e-Ts, (53) 
except that e and a, the step length, are now expressed in terms of m{T). With e = O ( 171 ^ ) 

AT 
m(T) 

T 2\ 



and s = 0(m(T)) the condition U = becomes = 0(1) p3"f , or 

0(A). (54) 



\ T 2 

\ c / 

This is the Ginzberg criterion for the onset of large fluctuations in the vicinity of T c , at 
which (51) becomes an unreliable basis for calculation. Thus, the Hagedorn transition of 
the previous section becomes subsumed in the conventional phase transition for the quantum 
field and the 'classical' and quantum pictures are reconciled. Our comments on attractive 
forces enforcing a maximum density if they are not to lead to unstabilities also survives. 
Moreover, although M 2 (T) diminishes towards the transition, so does the scalar dual field 
coupling strength g 2 , which we assume to be 0(e) = 0(m 2 (T)). Thus the maximum density 
does not diminish in the way we might have expected from the classical picture. 



Nonetheless, given that the string density is so high in the naive Kibble mechanism |L9 
it is always possible that this bound is breached. For type I strings another possibility exists. 
A first order transition is characterised by the production of latent heat. In the context of 
dual string theory it will be manifest as a first order transition in the string density p. A 
naive model for dislocation melting in solids has been based upon action Sn, treated as a 
'quantum' theory (h — 1) p5| . Since melting is, by definition, first-order, the 'quantum' 
fluctuations of the ^4-field in (48) can indeed induce a first-order transition in the order 
parameter p if the coupling f 2 of the orientable forces is strong enough The inclusion 



8 There is no contradiction in Type-II strings having a first order transition once the strength of 
the string forces can be chosen arbitrarily, and are not fixed by the underlying field theory. 
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of non-orientable forces clouds the issue. Such a transition is extremely difficult to predict, 
requiring an analysis of the dual theory beyond the mean-field approximation hinted at 
earlier. Calculations have been performed jn|] that extend the work of |18| to include the 
temperature dependence of the mass m(T) and the coupling strength / for strings (extreme 
type II) with short-range orientable forces. They show that, should a first-order transition 
occur, it too will be buried in the Ginzberg temperature range. The effect of this is that 
the string density can be very much lower (e.g. a factor 10~ 3 ) than we would have had 
otherwise. Such a low density is unlikely to have any problems for stability. 

As a final comment, it might be argued that the use of static networks, even initially, is 
inappropriate. However, the difference between an ensemble of static strings and an ensemble 
of relativistic strings, while present, is not as much as might be thought. For example, 
relativistic Nambu-Goto loops are described in terms of right-moving and left-moving modes. 
This doubling of modes changes the power behaviour of the prefactor to the Hagedorn 
exponential in the counting of states, in comparison to a static ensemble?]. However, this 
change is exactly cancelled by integrating over the loop centre-of-mass momenta, which 
effectively makes their centres of mass static. In the same sense, the time-independent 
action (51), if used as a basis for string saddlepoints, would give rise to zero-frequency 
modes that could be interpreted as centre-of-mass loop dynamics. 

Of course, we have been very simplistic in our neglect of quantum fluctuations about the 
strings. This has been considered elsewhere (last reference, and does not derail our 
conclusions to date. However, it does compound our inability to be quantitative with the 
dual theory. 



'It is this power which determines critical indices. 
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VII. CONCLUSIONS 



We have shown how the interaction energy for U(l) cosmic strings depends on the pa- 
rameter space of the model as well as on the relative orientations of two interacting strings. 

One of the crucial assumptions in all simulations of the evolution of networks of cosmic 
strings is that they intercommute in all cases. This property, however, should depend both 
on the topology of the problem which ensures the exchange of ends of the two strings 
and is model independent, but also on the dynamics of the system under the effect of its 
interactions. 

Our analysis reveals that type I strings experience attractive and, in the limit where 
the contribution from the gauge sector is negligible, non-orientable interactions. This is 
corroborated by several numerical studies and implies that in the absence of a strong effect 
at the superposition of their scalar cores two such strings will form a bound state of higher 
winding number, at a sufficiently low energy collision and for small relative angles. Such 
higher winding number states are well-known to be stable. A more detailed study concerning 



this problem is presented elsewhere [^6J , where we attempt to quantify the circumstances in 
which such configurations may occur. 

Type I strings are especially interesting for scenarios of thermal string production after 
a period of inflation, since the critical temperature associated with the phase transition is 



lowered by taking b >> 1 ||28|| . If the intercommuting of type I strings could be shown 
to be very inefficient such scenarios would have to be ruled out on the grounds of being 
cosmologically unacceptable. The most likely effect of the interactions, however, would 
probably be that of modifying the evolution parameters of a network of strings to some 



extent without jeopardizing the approach to a scaling regime pq| . 

Type II, global and critically coupled strings in our picture would intercommute. This 
is absolutely consistent with several numerical studies |||27 . 



Finally the knowledge of explicit forms for the interaction energy of strings is a crucial el- 
ement for the construction of a realistic statistical description of strings, as well as of vortices 
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away from critical coupling. In the last sections we have indicated, through the methods 
of a dual field theory, how strings with attractive forces do not destabilize the initial string 
network once the density becomes low enough. Observed originally for classical strings, it 
follows equally for strings produced by equilibrium fluctuations. Further, we see why b = 1 
ceases to be critical once arbitrary string orientations are taken into account. Of course, 
strings produced by fluctuations have to be frozen in by subsequent out-of-equilibrium de- 
velopment, and this has been omitted from our discussion. The whole mechanism of string 
production from a quantum theory, both with and without initial approximate thermal equi- 
librium, is considered elsewhere ||29|| . Nonetheless, the mechanisms proposed by Kibble |J 
for string production from fluctuations, from which we have been quoting, seem substantially 
correct. 
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FIGURE CAPTIONS 

Figure 1 : A schematic example of a two-string configuration and of the relevant quan- 
tities necessary to compute the corresponding interaction energy. 

Figure 2 : The field as measured by an observer away from a circular string loop. 
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